Controlled phase gate for solid-state charge-qubit architectures 
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We describe a mechanism for realizing a controlled phase gate for solid-state charge qubits. By 
augmenting the positionally defined qubit with an auxiliary state, and changing the charge distri- 
bution in the three-dot system, we are able to effectively switch the Coulombic interaction, effecting 
an entangling gate. We consider two architectures, and numerically investigate their robustness to 
gate noise. 
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I. INTRODUCTION 

The search for a workable quantum information pro- 
cessor is an effort that has captivated the attention of 
researchers in many disciplines. A quantum computer re- 
quires individual quantum logic elements, usually qubits, 
and entangling interactions between these elements 0,0. 

Solid-state proposals are widely seen as being some 
of the most attractive from the point of view of con- 
structional scalability, i.e., the ability to replicate many 
qubits. Furthermore, schemes compatible with present 
semiconductor technologies 0, 0, IE Hi are especially at- 
tractive because of their potential to leverage the associ- 
ated industrial experience |7], |8| • 

In this paper we concentrate on charge-based archi- 
tectures. Such systems were amongst the first pro- 
posed for quantum computing [fl and numerous versions 
have evolved recently ja, 11(1 ll ll . We are attracted 
to charge-based systems for three reasons: (1) proven 
high-fidelity readout compatible with single-shot opera- 
tions [l3|h (2) potential for high-speed (~ picosecond) op- 
erations |(| ; and (3) the ability to define variable dimen- 
sionality Hilbcrt spaces by appropriate partitioning |14| . 
The usual mechanism for coupling charge qubits is via the 
Coulomb interaction. In general, this interaction cannot 
be controlled without a variation in the charge distribu- 
tion in the qubits. In this paper we specifically address 
this issue and show how to make a scalable controlled 
phase gate that makes use of the Coulomb interaction. 

The Coulomb interaction is insensitive to minor vari- 
ations in the distance between the quantum dots and 
its strength lends itself to high-speed entangling oper- 
ations; conversely, its long-range nature makes it diffi- 
cult to effectively modulate interactions between qubits. 
Most charge qubit schemes so far proposed implicitly 
rely on a fixed, always-on Coulomb interaction between 
qubits 0, [T(J, HH Such schemes would necessarily require 



global control techniques 0,0,0 , which may be prob- 
lematic given the strength of the Coulomb interaction. 

In the following we first discuss local qubit operations 
(Sec. [nj an( i introduce the two-qubit interaction Hamil- 
tonian (SecHTHl- Tncn 

we describe how to realize a con- 
trolled phase gate and analyze its operation in terms of 
dynamic and geometric phases (Sec. IIV(I . Finally, we 
discuss various practical issues such as the detrimental 
effects of finite rise and decay times on the gate fidelity 
and how to correct them (Sec.^l, gate implementation in 
the presence of practical constraints on pulse lengths and 
tunnelling rates (Sec. lVI|) . and the effect of noisy controls 
(Sec. lVIljl and imperfect architectures (Sec lVIlH) on the 
gate performance. 



II. SINGLE QUBIT SYSTEM 

As will be shown below, a three-dot system is required 
to modulate the Coulomb interaction. We therefore sup- 
plement the canonical charge qubit with an auxiliary 
state and define a quantum element with three quantum 
dots and one charge, where the position of the charge on 
two sites defines the qubit and the third site defines an 
auxiliary state used for two-qubit interactions. We work 
only with the ground state of each dot. We further as- 
sume that the ground states are energetically close and 
sufficiently separated from higher-lying excited states, so 
that the excitation of these states can be neglected. The 
system can thus be approximated as a three level system 
with Hamiltonian 



H = J2^\d)(d\ 
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where \d) denotes the ground state of the electron local- 
ized in dot d for d — 1, 2, 3, ta the energy of state \d), Udd' 
the tunnelling rate between dots d and d' (for d' ^ d), 
and Xdd' = \d) (d'\ + \d') (d\. We shall assume that we can 
control the ground state energies e ( i and the tunnelling 
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For example, a Hadamard transform on the qubit sub- 
space corresponds to 
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It is possible to optimize the implementation of some lo- 
cal unitary operations by simultaneously changing mul- 
tiple control parameters. See appendix 1X1 



FIG. 1: Schematic of a single qubit-plus-auxiliary unit (bot- 
tom) and electronic potential (top). The large circles repre- 
sent the quantum dots, the black dot the shared electron. S 
and B surface electrodes serve to shift the ground state energy 
of the dots and change the height of the tunnelling barriers 
and thus the tunnelling rates between them. 

rates fidd't e.g., by varying the voltages applied to several 
control electrodes as illustrated in Fig. 

To implement local unitary operations, we inhibit tun- 
nelling to the auxiliary site by raising the barrier be- 
tween dots 2 and 3 (and/or 1 and 3 if applicable), or 
increasing the ground state energy of the auxiliary dot 
3. In practice, the precise functional dependence of td 
and Udd' on the control voltages applied should be deter- 
mined experimentally using Hamiltonian identification 
techniques |l9|. When //13 = /123 = the local Hamilto- 
nian can be rewritten as 

H = £l |l)(l| + e 2 |2)(2| + h 4 i 12 X 12 + e 3 |3)(3|, (2) 

and we can realize arbitrary unitary operations on the 
qubit subspace by changing the energy levels Cd and the 
tunnelling rate [i\ 2 . For instance, shifting the energy level 
of dot 2 by e 2 (i) for to < t < t\ results in a local phase 
rotation 

l> 2 (0) = |l)(l| + e-^|2)(2| + |3)(3|, (3) 

with cj) — €2(t)/Hdt. Effecting a tunnelling rate of fii 2 
for t < t < fi gives 

U a (a) = cos(a)fi2 - ism(a)X 12 + |3)(3|, (4) 

where the rotation angle is a — ji^ (t) dt and i\ 2 = 
|1)(1| + |2)(2|. By combining two phase rotations on the 
2nd dot with a tunnelling gate between sites 1 and 2, for 
example, we can implement any local unitary operation 
on the qubit subspace modulo global phases [3 

U(fa,a,fa) = U 2 (0 2 )U 12 (a) U 2 (fa), 



III. INTERACTION HAMILTONIAN 

To achieve entangling operations it is necessary to 
change the charge distribution of adjacent elements. Al- 
though almost any interaction will lead to entangle- 
ment 20], the resulting dynamics may be neither easy 
to utilize, nor robust against noise. 

Rather than choose an arbitrary interaction, it is more 
useful to consider geometries where the action of the 
Coulomb force on the qubit states is trivial. One way 
this might be achieved is by using shielding to eliminate 
direct interactions between qubits but not between the 
auxiliary sites as shown in Fig. [3 (a). Such a system could 
be fabricated in a GaAs 2DEG system fill l2l| or a pil- 
lar system 22], for instance. An alternative is to choose 
a geometry for which the effect of the Coulomb inter- 
action on the dynamics of the qubit subspace is trivial, 
resulting only in a global phase factor, as in Fig. [21(b). 
To fabricate such a structure, we would propose an ex- 
tension of the atomic placement techniques used for Si:P 
systems [23|, |24| . Both designs have the advantage of al- 
lowing the implementation of controlled two-qubit gates 
using simple pulse sequences. Note that using auxiliary 
sites has also been proposed as a way to allow generalized 
readout of quantum information in different bases |25| . 

The Hilbert space for two qubit-plus-auxiliary units 
is spanned by the states \dd') for d, d! = 1,2,3, where 
\d)k denotes the dth basis state for the fcth unit, and 
\dd') = ® \d') 2 are the tensor product states as usual. 
The total Hamiltonian of this system is 

H = H {1) ®h + h®H (2) +H C , (7) 

where ^3 is the identity matrix in dimension three, 
for k = 1 , 2 is the local Hamiltonian for the fcth qubit- 
plus-auxiliary unit specified in Eq. l[T|l. and He is the 
Coulomb Hamiltonian 

3 

Hc= 7dd'\dd')(dd'\. (8) 

d,d'=l 
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grounded metal barrier 




FIG. 2: (a) Two-qubit system with two auxiliary sites comprised of six quantum dots (filled circles). A grounded metal barrier 
shields the Coulomb interaction between nearby quantum dots except between the two auxiliary sites (3). The thick (red and 
blue) lines indicate the surface control electrodes and Sj , respectively, (b) 3D geometry for a two-qubit system with 

auxiliary sites, for which the Coulomb interaction between the two charges is constant if both are confined to their respective 
qubit subspace comprised of dots 1 and 2, but differs from the Coulomb coupling between the auxiliary dots provided that 
V46 2 + 2c 2 / 2a. Note that the inter-dot distances satisfy Hs^ - s^|| = V4& 2 + 2c 2 and Hs^ 1 ' - s^ 2) || = Hs^ 1 ' - s!j 2) || = 
^{a + b) 2 + c 2 for d, d! = 1, 2. 



The Coulomb interaction energies "fdd' are given by 

i^=^-yp- s f\\-\ o) 

where denotes the position of the dth quantum dot 
in the kth qubit-plus-auxiliary unit, e is the applicable 
dielectric constant, and e is the electron charge. In pure 
silicon we have e = 11.8eo, where eo is the dielectric con- 
stant in vacuum. 

If two sites and s# are separated by a sufficiently 
thick, grounded metal barrier then j^d' = 0. Hence, the 
Coulomb interactions for the shielded 2D geometry in 
Fig. [3(a) effectively vanish except for 733. Similarly, for 
the 3D geometry shown in Fig. [21(b), symmetry implies 

711 = 712 = 721 = 722 = 71 1 

713 = 723 = 731 = 732 = 72- 

We can cancel the effect of the Coulomb interaction be- 
tween qubit states by applying suitable bias voltages to 

(k) 

the energy shift gates to offset the energy levels e\ and 

4* b y -7i/ 2 and 4 fc) b y 7i/ 2 - 72 fo r k = 1,2, Thus, 
the Coulomb interaction Hamiltonian becomes 

ffc = 7eff|33)(33|, (fO) 

where 7 e g is the effective Coulomb coupling between the 
auxiliary states, i.e., He acts trivially on the system ex- 
cept if both electrons are in the auxiliary states. 

For a 2D architecture with shielding 7 c ff will usually 
be less than the free-space Coulomb interaction 733 due 
to screening and image charges. For a 3D geometry with 
energy offsets 7 e ff will be less than or equal to 733 — 272 + 
71 , with equality if the screening due to control electrodes 
etc. is negligible. 



For convenience, we choose the effective Coulomb en- 
ergy 7 e ff between the auxiliary sites as the unit of energy 
such that the free evolution Hamiltonian of the two-qubit 
plus auxiliary system is Hq = |33)(33|. All tunnelling 
rates are given in units of j e s/H and the canonical time 
unit is h/j e ff. 



IV. CONTROLLED TWO-QUBIT PHASE GATE 

To motivate the design of a two-qubit phase gate, we 
note that the phase acquired by an electron tunnelling be- 
tween two quantum dots depends on the energy difference 
between them. In our system the energy differences are 
determined by a combination of energy bias voltages and 
the Coulombic interaction between the auxiliary sites. 
Hence, by adjusting the voltage on the energy bias gate 
S% , for instance, we can shift of the energy of state |2)i 
such that the energy differences between states |21) and 
|31), and 1 23) and |33), respectively, are equal in magni- 
tude but opposite in sign, as illustrated in Fig. [21 In this 
case the phase acquired by a charge tunnelling between 
dots 2 and 3 of qubit I will have the same magnitude 
regardless of the population of state |3) 2 , but the latter 
will determine its sign. This observation is central to the 
operation of our controlled phase gate. 

To achieve a maximally entangling gate, the acquired 
phases in both cases must differ by an integer multiple of 
7r. Finally, except for the acquired phase, the charge must 
return to its initial state |2)i at the end of the tunnelling 
process, i.e., no population must remain in the auxiliary 
state 1 3)i. Combining all of these requirements leads to 
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FIG. 3: Energy level configuration during the second step: 
By applying an energy bias of e^ 1 ' = 1/2 we ensure that the 
energy gap between states 1 23} and |33), and |21) and 1 31) is 
±1/2, respectively. 



the following expressions for the tunnelling rate 



^23 



2n 



2k -1 



1, 



(11) 



for suitable positive integers n and k satisfying 2n > 
2k — 1, and tunnelling time 



47m 



T-2 



16 



2vr(2fc- 1). 



(12) 



1 



A detailed explanation of these results is provided in ap- 
pendices El and [D] Furthermore, the phase acquired by 
state \2l) is [2n - (2k - l)]vr/2, while that of state |23) 
is [2n + (2k — l)]7r/2, and hence the phase difference is 
(2k — 1)tt = 7T as desired. 

In the absence of constraints on the tunnelling rates 
and pulse lengths, the gate operation time is optimized if 
we choose n = k = 1 and ^3 = \/3/4. For details about 
how to choose n and k when there are constraints on the 
tunnelling rates and pulse lengths, see appendix lEl 

More explicitly, a maximally entangling controlled two- 
qubit phase gate 



f/phasc = |11)(11| + |12)(12| 

= diag(7 2 ,ii2)- 
can be realized as follows: 



|21)(21|-|22)(22| 



(13) 



1. Acting on the second qubit, swap the populations 
of the states 1 2)2 and 1 3)2 by lowering the tunnelling 

barrier between the 2nd and 3rd quantum dot for 

(2) (2) 
time Ti = 7r/(2/i2 3 ), where /iJ, 3 is the tunnelling 

rate. 



2. Acting on the first qubit, simultaneously raise the 

ground state energy of the 2nd dot by e^ 1 =1/2 
and lower the tunnelling barrier between the 2nd 
and 3rd dot to achieve the tunnelling rate given by 
Eq. lfTT|) for the time specified by Eq (jT2|l . 
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FIG. 4: Schematic representation of the phase gate operation 
for three initial basis states (a) |12), (b) |21), and (c) |22). 
The operation for |11) has been omitted because it is trivial. 
The filled (red) dots show the positions of the electrons at 
the beginning of the first step (top), during the second step 
(middle), and after completion of the third step (bottom). For 
initial configurations 1 21} (b) and 1 22} (c), the first charge is in 
a superposition of states |2}i and |3}i during the second step, 
and hence acquires a phase conditional on the population of 
state 1 3) 2, which is indicated by by lighter (pink) shading. 
Notice, however, that the electron starts and ends in state 
|2}i, except for the conditional phase acquired. 



3. Acting on the second qubit, repeat the first step to 
swap the populations of states 1 2)2 and 1 3)2 back. 

4. Acting simultaneously on both qubits, shift the en- 



ergy of states 1 2)i and |2) 2 by e 



(i) _ 



?7t/t4 and 

77/ V4, respectively, for some time 74, where 
£ = 1/2 - (n + k) mod 2. (See appendix|D]) 

The first three steps are illustrated in Fig. and Fig. 
shows results for a simulation of the gate for ideal pulses 
and no constraints. 

Since this gate combined with local unitary operations 
as described in Sec. [H] is universal, we can implement 
any desired two-qubit gate. For instance, a controlled- 
NOT gate is performed simply by applying a Hadamard 
transformation on the second qubit before and after the 
pulse sequence above. 

Further analysis shows that the conditional phase ac- 
quired by state |2)i is exactly twice the geometric phase 
acquired by the two-level subsystem Si = {|21),|31)} 
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2 4 6 8 10 12 14 16 

time [units of h/(2)«Y J] 

FIG. 5: Operation of the proposed controlled phase gate for 
ideal controls with no constraints. Shown are the control pa- 
rameter settings (a) and the corresponding evolution of the 
initial states (b) |12), (c) |21) and (d) |22). The evolution of 
1 11), being trivial, has been omitted. In each case, only the 
relevant populations and acquired phases are plotted. 



and S2 = {|23), |33)}, respectively. To see this, note that 
the Hamiltonians for S± and S2 are 



H Sl = j/2 + ^cr x 
3 - 

H S2 = -^h+ [i-o x 



1 



(14) 
(15) 



with fi = ( y/ (2n) 2 /(2k — l) 2 — 1) / 4, and we can visualize 
their evolution on the Bloch sphere as in Fig. Let 
s = (s x , s y , s z ) with s a = Tr(a a p) for a £ {x, y, z} be the 
usual Bloch vector. The initial states |21) and 1 23) for 
Si and 52, respectively, correspond to the Bloch vector 
so = (0,0,1), and their evolution in R 3 to a rotation 
about the axes di = 2(^,0,1/4) and d 2 = 2(^,0,-1/4), 
respectively. 

The^ pure-state non-adiabatic, cyclic geometric 

phase [2g is ^ 6Com = ^ = 27m (1 — cos 9), where fl is the 
solid angle subtended by the Bloch vector, n is the num- 
ber of times the Bloch vector rotates around the axis d, 
and 9 is the angle between the initial state vector So and 
the rotation axis d. Both states rotate with the same 



— Path of |21> 

— Path of |23> 




FIG. 6: Trajectory of the Bloch vectors associated with the 
two-level subsystems Si and S2 for ideal phase gate. For 
n = k = 1 the Bloch vectors si(t) (cyan) and S2(t) (black) 
rotate about di and cb, respectively, simultaneously complet- 
ing a single closed loop on the surface of the Bloch sphere and 
sweeping out solid angles of it (area "inside" the blue loop) 
and 3-7T (area "outside" the black loop), respectively. Hence, 
the difference between the areas is 2-7r, and the conditional 
phase acquired by state |2)i depending on whether state 1 3) 2 
is occupied will differ by n. 



frequency ||d|| = ||d'|| = (y/l + 16/i 2 )/2 = n/(2k- 1) 
and hence the relative geometric phase depends on 6 
Noting that 



cos 







di • Sq 2k — 1 



|di| 



2n 



(16) 



we see that the respective geometric phases acquired by 
Si and S 2 are 27m [1 - {2k - l)/(2n)] = n[2n - (2k - 1)] 
and 27m[l + (2fc-l)/(2n)] = 7r[2n+ (2k- 1)], i.e., exactly 
twice the conditional phases acquired by states |21) and 
|23). 



V. REALISTIC CONTROLS AND 
CORRECTION OF SYSTEMATIC ERRORS 

So far we have assumed ideal, piecewise-constant con- 
trol Hamiltonians. In reality, however, the pulses will 
have finite rise and fall times, for example, as shown 
in Fig. A naive implementation of the above scheme 
would perforce introduce systematic errors. These errors 
are evident in Fig. 0(1), which show the evolution of the 
system for the same parameter values as in Fig.[S]except 
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FIG. 7: Realistic square- wave pulse of amplitude A and length 
r with rise and decay time t s , modelled using error functions. 
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TABLE I: Gate error £ — 1 — T as a function of the pulse 
rise and decay time r s for non-ideal, uncorrected (£ u ) and 
corrected pulses (£ c ) for simulations with time step At — 
0.005. 



for pulse rise and decay times of r s = 1 time unit. Com- 
paring the trajectories shows that there are significant 
population and phase errors. In particular, the popula- 
tions of the states |13), 1 31} and 1 33) do not return to 0, 
i.e., population is lost to the auxiliary states, a poten- 
tially difficult error to correct. 

To quantitatively compare the gates implemented to 
the ideal gate we consider the average gate fidelity |27U28| 



J 7 =<* in |l/tpoutl7|*i 1 



(17) 



which is a measure of the overlap of the final state p ou t 
with the desired target state U^in) averaged over all in- 
put states. Since /5 ou t may extend to the auxiliary quan- 
tum dots, the fidelity includes the effect of population 
losses to these states. 

In the example shown in Fig. [5] pulse rise and decay 
times of one time unit reduce the average gate fidelity 
from one in the ideal case to ~ 73 %. Table [I] shows that 
even relatively small rise and decay times of the fields 
tend to reduce the average gate fidelity noticeably. This 
increase in the gate error appears to be mainly due to 
population loss to the auxiliary states. 

The first step toward improving the results is to re- 
alize that rise and decay times reduce the total pulse 
area. An ideal square-wave pulse (r s = 0) of duration r 
with amplitude A has a pulse area of At, while that of 
a similar pulse with rise and decay time t s < r/2 is only 
A(t — t s ) |2|j. Since the pulse area is an important con- 
trol parameter, we must therefore adjust either the field 
strengths or the pulse lengths, or possibly both. For the 
1st, 3rd and 4th step of the CPHASE gate correcting the 



pulse area is all that is needed to eliminate systematic er- 
rors due to non-zero pulse rise and decay times, and we 
can achieve this either by increasing the field strength or 
the pulse durations. 

During the crucial 2nd step, the energy of state |2)i 
should be kept constant = 1/2 while tunnelling is 

enabled, and the value of the tunnelling rate /Zgjj should 
be close to one of the desired values. Increasing the field 
strengths to compensate for rise and decay times is thus 
not an option for this step. For simplicity, we shall there- 
fore correct the pulse areas by increasing the duration of 
each pulse by r s , i.e., setting r' k = + t s , where Tu is 
the duration of the corresponding ideal pulse. 

To improve the accuracy of the 2nd step further, we 
raise the energy at least r s time units before we en- 
able tunnelling between states |2)i and |3)i, and lower 
it only after the tunnelling barrier has been raised again 
to inhibit tunnelling. Since changing the energy of state 
|2)i is a local operation on the first qubit, it commutes 
with the swap operation on states 1 2) 2 and 1 3) 2 - We can 
therefore begin to raise the energy of state |2)i before 
the 1st swap operation is completed, and we can begin 
the 3rd swap operation before £3 (t) has reached zero. 
However, the increased duration of the energy shift will 
induce an additional local phase shift, which has to be 
taken into account in the final step to achieve the desired 
gate. Let t\ = t[ — t s and t% = t[ + t' 2 + t s and 



-2 1} (*) dt- -T 2) 



(18) 



where Ti is the time required to complete the second swap 
operation for ideal, piecewise constant control pulses. 
Then the local phase rotation on the first qubit required 
in the final step is $ = lit — A<f> instead of in, where 
I = I — (n + k) mod 2 as before. Thus, in the final step 
the energy of state 1 2) x must be raised by 

(1) _ Til - A(j) 



(19) 



where 74 is the duration of the step for ideal, piecewise 
constant controls. 

Numerical simulations indicate that these corrections 
greatly enhance the performance of the phase gate. For 
instance, Table [I] shows that the average gate error £ c 
with these corrections is less than 10 -4 for all rise and 
decay times r s [3j| and Fig. [8] (II) shows that the trajec- 
tories very closely match those of the ideal gate. 



VI. GATE OPERATION TIMES AND 
PHYSICAL CONSTRAINTS 

In the hypothetical case of no constraints on the tun- 
nelling rates, energy shifts and pulse lengths considered 
so far, the first and third swap operation as well as the 
local phase rotations could, in principle, be implemented 
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FIG. 8: Control parameter settings (a) and evolution of the initial states (b) |12), (c) |21) and (d) |22), as well as (e) evolution 
of the two-level subsystems Si and 52 on the Bloch sphere under the operation of the controlled phase gate (n = k = 1) for 
controls with rise and decay time t s = 1 [ft/7 c ff] without correction of systematic errors (I) and with corrections described (II). 
Due to the rise and decay time of the controls, the trajectories of the Bloch vectors in (I) do not form closed loops, and the 
system does not return to its initial state (modulo the desired phase change). However, when the rise and decay times of the 
pulses are taken into account, population loss to the auxiliary states is negligible and the trajectories of the Bloch vectors again 
form closed loops (II). 



8 



arbitrarily fast and the gate operation time would be lim- 
ited mainly by the time required to complete the second 
step, i.e., r « T2 — 2n time units. In reality, however, 
the gate operation time is usually limited by physical and 
technical constraints that impose, for instance, a min- 
imum pulse length T m ; n and maximum tunnelling rate 
Mmax between the qubit and auxiliary quantum dot. 

To explore the consequences of finite tunnelling rates 
and switching speeds, let us consider a specific example 
of a 2D charge qubit architecture as shown in Fig. [2] (a) 
with auxiliary quantum dots spaced about 2a — 170 run 
apart j34|. Neglecting screening effects, the Coulomb en- 
ergy in this case is 7 c fr ~ 0.718 meV, and the charac- 
teristic time scale fry^ 1 ~ 1 ps 35] Thus, theoretically, 
two-qubit gate operation times of less then 20 ps could 
be achieved, as the previous sections show, and if one 
assumes that local operations can be realized arbitrarily 
fast then gate operation times of less than 7 ps would be 
possible. 

However, if the pulse lengths must be > 50 ps - 
about the limit of what is accessible with current tech- 
nology — then Eq. (|E1(I shows that the best possible 
choice of the parameters n and k that satisfies this 
constraint while minimizing the gate operation time is 
k = n = 5. The total gate operation time in this 
case increases to 3 x 50 + 187r 206.55 ps — as- 
suming we can achieve the required tunnelling rate of 

M23 = (l/4)v/(10/9) 2 - l[7cff/ft] ~ 1.21 x 10 11 rad s" 1 . 
The evolution of the system for the resulting gate is 
shown in Fig. (I). 

If the maximum tunnelling rate that can be achieved 
(e.g., without applying control voltages that might lead 
to a breakdown of the oxide layer separating the silicon 
substrate and the surface control electrodes, or pulses 
that might result in population loss to higher-lying ex- 
cited states) are lower than this, then the gate operation 
time will be increased further. For example, if we must 
satisfy [j$ < 10" 11 rad s' 1 then Eqs. l|E4| . (|E5)l show 
that the best choice for n and k is n = k = 7, which 
yields /i^j ~ 0.99 x 10 11 rad s -1 , and the gate operation 
time increases to 231.68 ps, as shown in Fig- El (H) - 



VII. GATE ROBUSTNESS FOR NOISY 
CONTROLS 

Another important issue in practice is the robustness 
of the gate for noisy controls. To study this effect we 
artificially added noise to our controls. Fig.llOlshows the 
increase of the gate error £ as a function of the magnitude 
of the noise added. The simulations were performed with 
a fixed time step of At = 0.005, for which the average 
gate fidelity in the absence of noise exceeded 0.9999. The 
noise functions r](t) were bounded 

\v(t)\<Vo (20) 



with Fourier transforms f]{oj) satisfying 

!0 u> = 

K/uj LU<LU (21) 
K/cj Ul > UJ . 

where A" is a constant. This type of noise corresponds to 
1// noise with a low-frequency cut-off, which is common 
in electronic devices. 

Although the addition of noise necessarily increase the 
gate error, the gate appears to be quite robust to this 
noise, as a typical example of the evolution of the system 
for noisy controls (Fig. Illfl shows. Not unsurprisingly, 
the curves for varying threshold frequencies ojq in Fig. 1101 
suggest that the gate is more sensitive to low-frequency 
than high-frequency noise, as the high-frequency compo- 
nents tend to cancel. 



VIII. GATE ROBUSTNESS FOR IMPERFECT 
GEOMETRIES 

A final issue that must be considered is the effect of 
manufacturing tolerances, which result in deviations of 
the qubit register from an ideal specification. In solid- 
state charge-based architectures of the types considered 
herein, the main issues appear to be imperfections in the 
placement or geometry of the quantum dots. The former 
is believed to be especially pronounced for quantum dots 
based on donor impurities in a Silicon matrix, while the 
latter will be relevant, e.g., for manufactured Ga/GaAs 
heterostructure quantum dots. As modelling all the pos- 
sible effects of imperfections in the fabrication process 
for heterostructure quantum dots would far exceed the 
scope of this paper, we shall focus on modelling the ef- 
fect of imperfect placement of the quantum dots on the 
gate performance. 

For solid-state architectures based on donor impurities 
(e.g., Phosphorus) in a Silicon matrix the accuracy of 
placement of the donors is limited. In a shielded 2D ar- 
chitecture inaccurate placement of the donors will mainly 
cause variations in the tunnelling rates as the structure 
of the silicon lattice introduces spatial oscillations in 
the donor wavefunctions, and hence tunnelling rates l30l 
analogous to those seen in exchange systems 0, l32j . 
There will also be minor changes in the strength of the 
Coulomb interaction between the auxiliary sites. How- 
ever, the actual Coulomb coupling strengths and tun- 
nelling rates as a function of the control voltages ap- 
plied can, in principle, be determined using Hamilto- 
nian identification techniques similar to those described 
in Ref. ^H), for instance, and the control scheme can 
then easily be adapted to achieve the desired gate for the 
actual system. 

The situation is different for 3D architectures [cf. Fig. [2 
(b)] that rely on symmetry to effectively cancel the 
Coulomb interactions between the qubit sites. Manufac- 
turing tolerances in this case will result in asymmetries of 
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FIG. 9: Control settings (a) and corresponding evolution of the initial states (b) |12), (c) |21), and (d) 1 22) , as well as evolution 
of the two-level subsystems S and S' on the Bloch sphere (e) under the operation of the controlled phase gate when (I) 
minimum pulse length constraints of T m m > 50 ps necessitate n — k — 5, and (II) simultaneous pulse length and tunnelling 
rate constraints r m i n > 50 ps and /tz 2 3 — 10 -11 rad s _1 necessitate n = k — 7. Note the decrease in the angle 81 (increase in #2) 
between the rotation axis di (cb) and the (positive) z-axis compared to the n = k = 1 case, and the resulting greater difference 
in the areas swept out by the Bloch vectors si (blue) and S2 (black) in a single loop: n/5 versus 197r/5 in (I); n/7 versus 277r/7 
in (II). Consequently, five and seven loops, respectively, are necessary to archive an area difference that is a multiple of 2n, and 
the time required to implement the phase shift increases to 187T and 267T, respectively. Also note the slight distortion of the 
loops near the north pole due to the dynamic change of the rotation axis during rise and decay periods of the pulses. Although 
present for all choices of n and k, it is more pronounced for larger n and k. 
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(0 = 25 (units of 2n/T ) 
[0 = 50 (units of 2u/T ) 
o> = 75 (units of 2it/T ) 
to = 1 00 (units of 2it/TJ 




0.05 0.06 0.07 
Magnitude of noise n 

FIG. 10: Average gate error £ as a function of the magni- 
tude r/o of the noise for various threshold frequencies uq. The 
solid curves show the mean of Sand the error bars indicate 
the range of £ over 10 simulations. The data were plotted 
horizontally offset to improve visual clarity. 



the geometry, which cannot be compensated easily, even 
if they could be identified precisely, and perforce reduce 
the gate fidelity. To estimate the effect of such errors 
we have performed computer simulations. For a 3D ar- 
rangement of six quantum dots as shown in Fig. [21(b), we 
randomly perturbed the positions of all six quantum dots 
by up to four lattice sites in the x- and y-directions, and 
±1 monolayer in the z-direction, assuming a lattice con- 
stant of 0.3 nm, values that appear feasible with current 
fabrication technology |23| . and numerically computed 
the average gate fidelity T for each perturbed system. 
Note that it was assumed here that the tunnelling rates 
can be kept steady despite the effect of the silicon lattice 
by adjusting the control voltages appropriately. 

The numerical results in Table HP suggest that the ro- 
bustness of our phase gate with respect to asymmetries 
depends significantly on the choice of interqubit spac- 
ings and the distance between the auxiliary sites. Con- 
cretely, the data suggest that the robustness of the gate 
with respect to random pertubations of the geometry 
is maximized by minimizing the distance between auxil- 
iary sites while maximizing the distances between qubits. 
There also appears to be a strong relation between the 
robustness with respect to asymmetries and the effec- 
tive Coulomb interaction between the auxiliary sites in 
the system, suggesting that maximizing the later quan- 
tity might be advantageous. However, it should be noted 
that stronger Coulomb coupling implies shorter gate op- 
eration times unless the gate parameters arc changed, 
which might affect the gate fidelity. A final computa- 
tion for a (target) geometry with a — 20, b — 100 and 
c = 10 showed that the average gate fidelity over 100 
random perturbations ranged from 0.9944 to 1.0000 with 




* |< 22|22>| £ 
0.5 * l<23|23>| 2 \ 

a |( 33|33)| 2 x 
_ 1 . - angle[<y<0),y(i)>l/iT , . . | . . 
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time [units of h/(2jry )] 



Path of |21> 
Path of |23> 




(e) 



FIG. 11: Control settings (a) and corresponding evolution of 
the initial states (b) |12), (c) |21) and (d) |22>, as well as 
evolution of the two-level subsystems Si and S2 on the Bloch 
sphere (e) under the operation of our phase gate for noisy 
controls (770 = 0.1, luq — 50). The noise has no noticable 
effect on the evolution of the populations. The relative phases 
exhibit some wiggles but they mostly seem to average out 
over the duration of the gate. The average gate error was 
£ = 0.708 x 10" 3 . 
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a \ b 


30 


50 


70 


90 


20 


0.964022 


0.996631 


0.998153 


0.998855 


40 


N/A 


0.909071 


0.992692 


0.998410 


60 


N/A 


N/A 


0.831827 


0.989618 


80 


N/A 


N/A 


N/A 


0.808333 



TABLE II: Mean of the average gate fidelity T for different 
geometries with fixed intra-qubit spacing of c = 10 nm. For 
each target geometry the mean of T was computed for 30 
randomly perturbed systems. The gate parameters for all 
simulations were n — k — 1 and the time steps were chosen 
such that the average gate fidelity of each unperturbed systems 
was > 0.99995. Notice that the mean of T increases sharply 
for decreasing values of a, and noticably for increasing values 
of b. 



a mean (standard deviation) of 0.9989 (0.0013). These 
results suggest that even 3D charge-qubit architectures 
without shielding could be designed to be rather robust 
with respect to fabrication errors. 



IX. CONCLUSIONS 

We have presented two scalable achitectures for solid- 
state charge qubits that permit controlled entangling op- 
erations between pairs of qubits. We add that these op- 
erations could also be performed in parallel for cluster 
state preparation, a concept that we will explore else- 
where. The key feature of both geometries is that inter- 
actions between qubits are mediated via auxiliary quan- 
tum dots, while direct interactions between qubits are 
suppressed either through the use of shielding, or a three- 
dimensional design to cancel interactions. Controlled en- 
tangling gates can therefore be implemented by switching 
the charge distribution, and hence the Coulomb interac- 
tion, between qubits using the auxiliary dots. Both sys- 
tems should be realizable using existing or near-future 
fabrication techniques. 

In particular, we have shown explicitly how to realize 
a controlled phase gate, i.e., a universal, maximally en- 
tangling two-qubit gate by a simple four-step procedure. 
The crucial step in the sequence is the controlled tun- 
nelling between a qubit state and an auxiliary dot, which 
induces a phase shift conditional on the occupation of an 
adjacent (auxiliary) quantum dot. The scheme is suffi- 
ciently flexible to accommodate practical constraints on 
both pulse lengths and tunnelling rates. Analysis of the 
gate operation shows that the controlled phase shift can 
be explained in terms of dynamic and geometric phases. 
Due to the strength of the Coulomb coupling the gate 
operation is surprisingly fast. In the absence of pulse 
length constraints gate operation times on the order of a 
few picoseconds would be theoretically possible, and even 
with currently realistic constraints on the pulse lengths 
and tunnelling rates, gate operation times around 200 ps 
should be attainable. 

Using computer simulations, we have also studied the 



effects of systematic errors on the gate performance. The 
simulations show that finite pulse rise and decay times 
tend to result in population loss to the auxiliary states, 
and can noticeably reduce the average gate fidelity. The 
average gate error increases sharply with increasing pulse 
rise and decay times. The ideal scheme, however, can 
easily be generalized to compensate for such system- 
atic errors. Simulations show that these corrections can 
greatly improve the average gate fidelity; the corrected 
pulse sequences always achieved average gate fidelities of 
> 99.99%, and in some cases the average gate fidelity 
increased from ~ 70% without corrections to > 99.99%. 
These results suggest that (experimental) characteriza- 
tion of the rise and decay times of the pulses is very 
important to achieve high gate fidelity. 

We have also performed simulations to assess the ro- 
bustness of the gate to noisy control pulses. The re- 
sults suggest that the gate is quite robust with regard 
to both bandwith-limited l/f noise and white noise. 
High-frequency noise tends to effectively cancel. Low- 
frequency fluctuations of the control pulses can reduce 
the gate fidelity. Surprisingly, however, the resulting av- 
erage gate errors are generally very small compared to the 
systematic errors. In most of our simulations, the average 
gate error increased only from < 10~ 4 to < 10 -3 , even 
for very noisy controls. Our simulations further suggest 
that even a 3D geometry without shielding, which relies 
mainly on symmetry to cancel the effect of the Coulomb 
interaction between qubits, can be made to be quite ro- 
bust to misalignment errors during fabrication if the pa- 
rameters of the geometry are chosen sufficiently carefully. 
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APPENDIX A: LOCAL OPERATIONS 

We observed in Sec.|U]that any local unitary operation 
can be realized by concatenating three elementary rota- 
tions. It is worth noting, however, that some local uni- 
tary operators can be implemented more efficiently if we 
allow simultaneous changes of more than one control pa- 
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rameter such as the energy difference Aei 2 = (£2 — ei)/2 
and the tunnelling rate /i 12 between states |1) and |2). 
To see this note that we can rewrite Eq. (J2J) as follows: 

ff(e<j,Mi2) = ei2ii2 + Aei2Zi2 + Vi2^i2 + e 3 |3)(3| (Al) 

where e 12 = (e x + e 2 )/2 and Z 12 = |2)(2| - |1)(1|. Thus, 
if we apply constant energy shifts and effect a fixed 
tunnelling rate /ii 2 between dots 1 and 2 (while all other 
tunnelling rates are kept zero) for < t' < t then we 
generate the unitary operator: 



1 a? .sin(at) 
cos(at)ii2 — 1 -— — — 



(Aei 2 Zi2 + h^ii 2 Xi 2 ) 



exp(-iiei 2 /ft) + exp(-iie 3 //i)|3)(3| 



(A2) 



where a = W (Aei2/fr) 2 + n\ 2 - Thus, we can realize a 
Hadamard transform, for instance, in a single step by 
setting Aei2 = —fr^u 7^ for time t = 7t/(2\/2mi2)- 



APPENDIX B: FIRST SWAP OPERATION 

The total Hamiltonian for the first swap operation is 
Hi = H Q + I 3 <g> i? (2) where 



^ 2 )=^ 2 2 3 ) (|2)(3| + |3)(2|) 



(Bl) 



and Hq = 1 33) (33 1 as before. Applying this Hamiltonian 
for time < t' < t gives rise to the unitary operator 
Ui{t) = exp[-(it/h)Hi\. If we set t = r x = 71/(2^) we 
obtain the block-diagonal matrix 



U(n) = diag(l, -iX, 1, -iX, 1, W) 

' 1 x 



(B2) 



where X = \ and W is a symmetric 2x2 matrix 

[lo) 

with 



-in/2 



W 22 = e-"^ 2 



,-»ti/2 



4';s 
Ui 



cos 



V 2 / ui Sm V 2 



and Mi = Y 1 + 4[yU 2 3^] 2 . J7i (ti) swaps the population 
of states |2)2 and 1 3) 2 provided that state 1 3) 1 is not 
occupied. 



APPENDIX C: SECOND SWAP OPERATION 

The total Hamiltonian for the second gate is H2 = 
H + ® I3 where 

i?( 1 ) = i|2)(2| + ft M2 1 3 ) (|2)(3| + |3>(2|) (CI) 



and Hq = 1 33) (33 1 as before. Applying this Hamiltonian 
for time < t' < t gives rise to the unitary operator 
U%(t) = exp[— {it/K)H2[, which has the general form 





(h 








&(t) = 





A 


B 




Vo 


B 


C 



(C2) 



where A = diag(a, a, a'), B — diag(6, 6, b') and C = 
diag(c, c, c'). We must choose the gate operation time 
t 2 such that B = 0, i.e., b = b' = 0. Since 

(t)#-(t) <-> 
— t)^(t) <«) 



with w 2 = 
4n7r for some integer n, or 

7"2 



1 + lGf/i^] 2 , this is equivalent to u 2 T 2 = 



Ann 



'l + ie^] 2 

This choice of gate operation time gives 

£Mt 2 ) = diag(l, 1, l,a,a,a' ,a, a, a') 
where we have 

a = (— 1)™ exp(— inir/u2), 
a' = (— 1)" exp(— i3n7r/u 2 ). 



(C5) 



(C6) 



(C7) 



APPENDIX D: LOCAL PHASE ROTATIONS 

Composing the three swap operations, noting 1/3^3) — 
Ui(n), gives 

U = Ui(n)U 2 (T2)Ui(n) 

= diag(l, — 1, —1, a, —a', —a, a, W diag(a, a')W) 

(Dl) 

whose projection onto the two-qubit subspace 
{|11),|12),|21),|22)} is U = diag(l,-l,a,-a') with a 
and a as in (|C7|) . which is not quite a controlled phase 
gate yet. However, if we multiply U by ® U^- 2 \ 
where 



f/W = diag(l,(-l)"exp(wrn/it 2 )) 
= diag(l,-l) 



(D2) 



which corresponds to simultaneous local phase rotations 
on state |2)i and |2) 2 , we obtain 

U {1) ® U {2) U = diag (1, 1, 1, exp(-i27rn/u 2 )) (D3) 
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where m 2 = \Jl + 16 2 and 



M) - 



2n 



M23 ~ 4 VV 2k - 1 



- 1 



(D4) 



for positive integers n and k satisfying 2n > 2k— 1. Then 
clearly 27m = (2k — l)7ra 2 , i.e., we have 



pulses usually cannot be arbitrarily short and there the 
tunnelling rates cannot be made arbitrarily large, and we 
can accommodate these constraints by choosing suitable 
values for the parameters n and k. 

If the mimimum pulse length is r 2 > T m i n but the 
tunnelling rates are unconstrained then we simply set 



k 



U {1) <g> U (2) U = diag(l, 1, 1, -1) = £7, 



phase- 



(D5) 



1 



in 



(El) 



To achieve a controlled phase gate on the qubit space we 
thus require 



where \x] is the smallest positive integer > x, to optimize 
the gate operation time while satisfying the pulse length 
constraint. 

If we must satisfy fi^3 — f^max but there is no con- 



cxp(-i27rn/u 2 ) = -1 4^ 2nn = (2k - 1)ttm 2 . (D6) straint on the length of the control pulses, we set 



Substituting u 2 = 2n/(2k - 1) yields (-l)" e W" 2 = 
e in(n+k-i/2) ^ which ghows that fj(i) = u^ } (£tt) with £ = 

1/2 — (n + k) mod 2 and — (tt), i.e., we must 



k = IKu^-l)- 1 ] /2\+l 
n = [(2k - l)u max /2\ 



(E2) 
(E3) 



(i) _ 



7r/r4 and e 2 = tt/t4 for some time t = T4 in 



set e 



the final step. 



APPENDIX E: TUNNELLING RATE AND 
PULSE LENGTH CONSTRAINTS 

In absence of constraints on the tunnelling rate /i 2 ^ it 
is easy to see that the optimal choice of the parameters n 
and k in the second step is k = n = 1 , which implies 112 = 
2 and (X23 = V / 3/4- In practice, however, the control 



where u max — \/^§[£; nax + 1 and [x\ is the largest (pos- 
itive) integer < x, to satisfy the constraint and optimize 
the gate operation time. 

Finally, if we must satisfy r 2 > T m in and fj,^ < Umax , 
we choose 

k = max{fci, rr min /47r+ 1/2]} (E4) 
n = [(2k-l)u rnax /2\ (E5) 



where k x - 

Umax = 



l)" 1 ] /2j + 1 as in and 



max + 1 as before. 
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